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ABSTRACT 

Incomplete carbon-oxygen detonation with reactions terminating after burn- 
ing of C^^ in the leading + C^^ reaction (C-detonation) may occur in the low 
density outer layers of white dwarfs exploding as Type la supernovae (SN la). 
Previous studies of carbon-oxygen detonation structure and stability at low den- 
sities were performed under the assumption that the velocity of a detonation 
wave derives from complete burning of carbon and oxygen to iron. In fact, 
at densities p < 10^ g/cm^ the detonation in SN la may release less than a 
half of the available nuclear energy. In this paper we study basic properties of 
such detonations. We find that the length of an unsupported steady-state C- 
detonation is ~ 30 — 100 times greater than previously estimated, and that the 
decreased energy has a drastic effect on the detonation stability. In contrast to 
complete detonations which are one-dimensionally stable, C-detonations may be 
one-dimensionally unstable and propagate by periodically re-igniting themselves 
via spontaneous burning. The re-ignition period at p < 10^ g/cw? is estimated 
to be greater than the time-scale of a SN la explosion. This suggests that prop- 
agation and quenching of C-detonations at these densities could be affected by 
the instability. Potential observational implications of this effect are discussed. 
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Introduction 



Type la Supernovae (SNIa) are thermonuclear explosions of carbon-oxygen white dwarfs 
(CO- WD) in binary stellar systems. Explosion models of SNIa currently discussed in the 
literature include explosions of Chandras ekhar-mass WD, su ch as delayed detonation, pulsat- 



ing d e layed detonations and its variants (iKhokhlov 
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Jordan et al. 
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20061: 



Jackson et al. 



2008 



Meakin et al.l 



Bravo et al 
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Gamezo et al.lboosj Ihivne et al.l 



2010), gra vitationally confined (jPlewa et al. 



2009[) and pulsating reve rse detonations 



sub-C handrasekh ar W D JWooslev et al.lll98ol:lNomot Jl982l: 



Bravo fc Garcia-Senzl 120091). explosions of 



Livne fc Glasner 1991: Livne fc Arnett 



19951 : iFink et al.l 120071 . |2010| ) and super-massive WD (jPfannes et al.l 120101 ). These models 
differ in their assumptions about initial conditions, ignition processes, whether the explosion 
involves subsonic deflagration or not, and other details, and they have a varying success in 
explaining basic observations of SN la (this subject is outside of the scope of the paper). 
A common feature of the models is that all of them involve, in one way or another, the 
detonation mode of burning. 

Detonation is a supersonic wave of burning in which reactions are triggered by a strong 
shock. The energy released by burning maintains the shock strength and prop agation of 
the detonation wave . A steady-state det onation is described by the ZND theory (IZeldovich 
1940l : Ivon NeumannI Il942l : iDoringI Il943l ) according to which the detonation wave consists 
of a leading planar shock followed by a one-dimensional reaction zone. A ZND detonation 
structure is usually unstable due to a strong positive feedback between hydrodynamical 
perturbations and energy release inside the reaction zone, and a real detonation propagates 
non-steadily with strong oscillations and t ransverse waves which form a multi-dimensiona l 
cellular structure of a detonation wave (e.g. iFickett &: Davislll979l : lLedll984j : IShepherdll2009l ). 



Explosive burning in SN la consists of three distinct stages - burning of carbon to O, 
A^e, Na, Mg and some Si (C-burning), subsequent burning of oxygen and formation of 
Si-group elements (0-burning), and flnally burning of silicon to Fe-group elements {Si- 
burning). Nuclear statistical equilibrium (NSE) in burned matter sets in at the end of Si- 
burning. Time-scales of C, O, and 5'i-burning, tc, to and tsi, differ by orders of magnitude, 
tc « to « tsi, and increase exponentially with decreasing temperature. Due to the 
existence of three distinct stages of burning, a CO detonation waves in SN la consists of C, 
O, and Si-burning layers of increasing thickness, xc << xo << xsi, following each other. 
All three thicknesses depend on density, p, and increase exponentially with decreasing p. 
Thickness of Si and 0-burning layers, xsi and xo, become greater than the characteristic 
scale of an exploding WD (2:^ 10^ cm) at p < 10^ g/c m^ and p < 10^ g/cm^, respectively 



e.g. lKhokhlovlll989l : iGamezo et al.lll999t ISharpd 119991 ). Equivalently, it could be said that 
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at these densities tsi and to become greater than the characteristic explosion time-scale of 
WD ~ 1 sec. An important observational property of SN la, as was early inferred from 
near-maximum spectra, is the presence of intermediat e mass elements such as Si, S, Ca 
Mq, Na and O in outer layers of a SN la envelope (e.g. |Pskovskiilll969l : [Branch et al.lll982 



19831 ). Presence of these elements can be explained as a consequence of incomplete burning 
at low densities p < 10'' g/cm^. 

Whether and how the instability and cellular structure of a detonation influences the 
explosion of a SN la has been a long-standing question. Time-dependent interplay of nuclear 
reactions and hydrodynamical motions in a detonation wave may affect explosive nucleosyn- 
thesis and quenching of burning. Incomplete C-burning may leave some carbon intact. This 
may be important for prediction of chemical composition and formation of light curves and 
spectra of SN la before and near maximum light. Due to the resolution limitations these 
effects were studied in ideahzed const ant- density conditions representing small parts of an 
exploding white dwarf. A one-dimensio nal longitu dinal instability of C O detonations was 
studied numerically in iKhokhlovl (119931 ) (see also iKoldoba et al.l (119941)'). Cellula r struc - 
ture of CO detonat io ns was obtain e d in numerical simulations of Isoisseau et al. ( 1996 ): 
Gamezo et al.l (119991 ) ; iTimmes et al.l (120001 ) . Calculations show that cellular structure asso- 
ciated with O and 5'i-burning might be potentially important and affect burning at densities 
P < 10'' g/cm^ but cellular structure related to C-burning appears to be very small at all 
densities of interest, say, at p > 10^ g/cm^. 

For low-density detonations the simulations mentioned above contain a serious defi- 
ciency. The reason is as follows. The simulations use background conditions (density, tem- 
perature, etc.) and the detonation velocity as input parameters. These parameters are used 
to establish initial and boundary conditions for a simulation. The simulation then gives 
information about scales and stability of a detonation wave. In the above simulations it was 
assumed that the input detonation velocity corresponds to a complete burning to NSE. For 
detonations at low density this leads to a contradiction: Si, and 0-burning at low densities 
may not occur and a detonation wave may consist of a shock wave followed either by C 
and 0-burning layers (0-detonation) or by a C-burning layer alone (C-detonation). Nuclear 
energy released in such detonation waves will be smaller than the NSE energy release. It 
is this reduced energy release which must be used for calculating the input detonation ve- 
locity. The assumption of NSE will lead to a higher detonation velocity, over-prediction of 
post-shock temperature and density, under-estimation of a detonation wave thickness, and 
may affect detonation stability and propagation. 



In this paper we re-visit the problem of detonation structure and stability taking into 
account the reduction of energy release in incomplete detonation waves. We try to answer the 
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foUowing questions: What are the basic properties and characteristic spatial and temporal 
scales of incomplete O and C-detonation waves at low densities? How unstable are these 
detonations? What is the time-scale of the instability and how does it compare to other 
scales of the problem? To what extent may the instability affect quenching of a detonation 
wave and the resulting nucleosynthesis? In this paper we study ZND structure and one- 
dimensional stability of incomplete detonations. Propagation of a detonation wave is also 
influenced by mult i- dimensional effects but a one-dimensional investigation should be the 
first step. We find in this paper that a C-detonation could be a factor of ^ 30 — 100 thicker 
than previously thought and that C-detonations could be highly unstable with respect to 
one- dimensional longitudinal perturbations. Multi-dimensional study of C-detonations will 
be presented in the next paper. 

The paper is organized as follows. Sect. IHpresent governing equations and input physics. 
Sect. [3] discusses properties of steady-state detonations with complete and partial energy 
release. Sect. H] describes numerical simulations and presents the analysis of one-dimensional 
stability of the detonations. The results and potential implications for SN la modeling and 
observations are discussed in Sect. El Numerical method and test simulations are described 
in Appendix \^ 



2. Formulation of the problem 

Propagation of a one- dimensional detonation wave is described by the time-dependent, 
compressible, reactive flow Euler equations of fluid dynamics. 



dp _ d{pu) 

P\ 

(2) 



dt dx 
dpu d{pu^ + P) 



dt dx 

dE_ _ d{u{E + P)) 

dt dx 

dp? d{puY) 



+ pq, (3) 



dt dx 

for mass density p, velocity u, energy density E = p(e + ^), where e is internal energy per 
unit mass, and the composition of reactants Y = {Yi, Yjy} defined as Yi = rii/pNa, where 
is the number of reactants, - their number densities, and Na is the Avogadro number. 
R = {Ri, Rn} are the corresponding net reaction rates which are functions of Y, p, and 
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temperature T. The energy generation rate due to nuclear reactions is 



N 



q = Q-R = Y^QiRi 



(5) 



where Q = {Qi, Qn} are the binding energies of nuclei. Numerical method of integration 
of equations (II]) - (H]) is described in Appendix |Al 

Nuclear kinetics is described by the a-network which h as been extensively used in pre 



vious studies of detonatio n stability and cellular structure (IKhokhlovl Il993l : iGamezo et al. 



1999 



Timmes et al.ll2000h . The network consists of = 13 nuclei He"^, C^^, O^^ A^e^", 
^^28^ ^32^ ^^36^ Q^io^ rp-ii^ ^^48^ ^^52^ ^-56^ rj.^^ network takes into account 

binary reactions between a-particles and heavier nuclei C^^ + a -H- O^^ + 7, O^^ + a ^ Ne^^, 
Fe^^ + a ^ Ni^^; C^^ + C^^ + O^^, and O^^ + O^^ and a triple-alpha reaction 
3a o C^^ + 7. Effective reaction rates of binary reactions involving a particles are calcu- 
lated as the sums of contributions of 7), (a, n)(?T,, 7), (0,7) reaction channels in- 
volving p, n, an d 7-photons. Forwa r d reaction rate s and partition funct i ons a re taken from 
compilations of iFowler et al.l ( 1l978l ) ; IWoosley et al.l (119781 ) ; iThielemannI (119931 ) . Reverse re- 
action rates are calculated from the principle of detailed balance. The network captures the 
multi-stage nature of explosive CO burning and correctly reproduces the time-scales of C, 
O, Si burning stages, and the onset of NSQE and NSE in nuclear matter. The equation of 
state includes contributions from ideal Fermi-Dirac electrons and positrons with arbitrary 
degeneracy and relativism, equilibrium Planck radiation, and ideal Boltzmann nuclei. 



3. Steady-state CO detonations 
3.1. ZND detonation in CO mixtures 



Steady-state detonation solutions provide basic temporal and spatial scales of a detona- 
tion wave are also used as initial conditions for time- dependent numerical simulations. In a 
reference frame moving with the detonation velocity, D, all variables inside a one-dimensional 
steady-state detonation wave are a function of the distance, x, from the leading shock and 
do not depend on t. We denote with subscripts 0, s, and d hydrodynamical variables in 
unburned cold matter, in unburned matter located immediately behind the leading shock, 
and in detonation products, respectively. Hydrodynamical variables inside the reaction zone 
are related to hydrodynamical variables in unburned matter through the Hugoniot relations, 

pu = poD, (6) 

P-Po = -{poDf-{V-Vo), (7) 
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e-eo + -(P + Po)-(V^-K)) + g = 0, 
where = 1/p is the specific volume, and 

q = Q.{Y,-Y) 



(8) 



(9) 



is the energy released in a fluid element by burning from to x. On a P-V diagram (Fig. [T]) 
physical conditions of matter inside the reaction zone are located on the Rayleigh-Mikhelson 
(R-M) line ([7]) with the slop tg a = —{poD)'^. For g = 0, ([HD gives the usual shock adiabat, 
S. Post-shock conditions are determined by the intersection of the R-M line with S. For 
Q = Qd, & gives the detonation adiabat. Conditions in detonation products correspond to 
the intersection of the same R-M line with the detonation adiabat. In case of a nuclear 
statistical equilibrium (NSE) Yd = Yd^NSEiTii, Pd) is a unique function of temperature and 
density and qa = Q-Yd^NSE = Id^NSEiTd, Pd) is a function of Td and pd as well. The detonation 
adiabat with q = qd,NSE is marked on Fig. [1] as NSE. The evolution of chemical species 
variables Y with x is described by a steady-state form of (jl]), 



dY 
dx 



u 



Rip,T,Y). 



(10) 



Conditions inside the reaction zone are determined by intersection of t h e R-M line with 
partial adiabats ([8]) with q^O determined by (ZND; iFickett fc Davi^ Jl979h l 



The detonation solution si - CJnse in Fig. [T] is the Chapman- Jouguet (CJ) detonation 
regime with the minimal possible detonation velocity, D = Dcj- An unsupported detonation 
propagates in this regime if g is a monotonic function of x. Overdriven detonations are 
possible in this case for all velocities D > Dcj- This is the case for all detonations at 
po < 2 X 10'' g/cjD? considered in the paper. For higher densities the energy release in 
0.5C -|- 0.5O detonations is not monotonic and an unsupported detonation propagates in a 
pathological regime with some velocity D* > Dqj- Overdriven detonations in this case are 
possible for D > D*. For a few simulations of high-density detonations presented in this 
paper we used as initial conditions slightly overdriven detonation solutions with D very close 
to D*. A detailed explanation of a pathological detonation can be found i 



Khokhlov ( 


1989) and 


Sharpe ( 


1999) 



3.2. CO detonations with complete energy release 

As mentioned in Sect. [T|, explosive burning of CO mixtures proceeds through three 
consecutive stages - burning of C^^ in the leading C^^ + C^^ reaction, burning of O^^ and the 
formation of Si-group elements, and Si-burning which creates Fe-peak elements. The state 
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of NSE is reached at the end of S'i-burning. As a resuh, the detonation waves in carbon- 
oxygen have a muhi-layered structure with C-burning, 0-burning, and Si-burning layers 
with increasing thickness, xq << xq << xsi, which follow each other. The CJ detonation 
structure for a 0.5C -|- 0.5O detonation at po = 3 x 10^ g/cm^ is shown in Fig. |2] as an 
example. The distinct C-burning, 0-burning and S'z-burning layers are clearly visible. C- 
burning releases qc — 0.35 MeV/nucleon of nuclear energy, at the end of the 0-burning the 
total released energy is qo — 0.65 MeV/nucleon, and the total energy released at the end of 
the reaction zone is qd,NSE — 0.8 MeV/nucleon. The detonation products at this density are 
mostly composed of Ni. 

In what follows we will characterize detonations in CO mixtures with the overdrive 
parameter 

\^C.J,NSE 

calculated with respect to the CJ detonation velocity of the complete detonation in which 
burning proceeds to NSE. Pathological detonation occurring at po > 2 x 10'' g/cm^ has an 
overdrive /* = {D* / Dcj^nse)'^ > 1- Calculations show that deviation of /* from one are 
relatively small. For po = 3 x 10'' g/cm^ /* ~ 1.01. For higher densities /* increases to 
/* ~ 1.07. Fig. |3] shows an overdriven detonation solution for po = 10^ g/cm^ with / = 1.08 
which is very close to the pathological regime at this density. Again, the distinct C*-burning, 
0-burning and S'z-burning layers can be clearly seen. Similar to the previous case, O-burning 
releases qc — 0.36 MeV/nucleon of nuclear energy but due to photo-dissociation of nuclei 
in 0-burning and Sz-burning layers qo is reduced to ^ 0.5 MeV/nucleon and the overall 
detonation energy to qd,NSE — 0.37 MeV/nucleon, resulting in the non-monotonic energy 
release. For detonations with non-monotone energy release we find that maximum of q 
occurs near the maximum of the concentration of Si. The energy release in O-burning layer 
is always monotonic so that O-detonation will always be not pathological. 

Fig. m gives xc, xq, and xsi for unsupported 0.5O + 0.5O detonations as a function 
of Po; Xc and xq are defined as half-reaction zone thickness where the concentrations of 
O^^ and O^^ diminish to one-half of their initial values, respectively; xsi, is defined as a 
thickness where Si decreases to one-half of its maximum value inside the reaction zone. The 
total thickness of the detonation wave, x^se, corresponds to 90 % of the Ni synthesized 
in the detonation wave. The horizontal dotted line in Fig. H] shows the characteristic size 
of an exploding WD, Rwd ~ 10*^ cm. For p < Pnse — 10'' g/cm^ xjysE > Rwd which 
means that burning in a detonation wave will not have time to reach the state of NSE. At 
P < Psi — 5 X 10^ g/cm^ the reaction zone of a CJ detonation will consists of O-burning layer 
followed by O-burning layer but Si-layer will be absent. At p < po — 10^ g/cm^ xq > Rwd 
and the reaction zone will consist of O-burning layer alone. 
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3.3. Carbon detonations with reduced energy release 

Due to a large separation between the reaction scales of C, O, and S'i-burning, an 
incomplete detonation may be considered as an isolated C-detonation (or an 0-detonation) 
wave with reduced energy release qc — 0.35 MeV/nucleon (or qo — 0.65 MeV/nucleon for O- 
detonation). At densities p < 10^ g/cm^ photo-dissociation of nuclei is insignificant and qc 
and qo are practically density-independent. For C-detonation the situation is schematically 
illustrated in Fig. [1] The figure shows a C-detonation adiabat with energy release qd = qc- 
The R-M line s3-CJc corresponds to a CJ C-detonation with velocity Dcj,c < Dcj,nse- 

By analogy with f llip we can also characterize C-detonations with the overdrive param- 
eter defined relative to Dcj,c, 

Then C-burning layer of a CJ NSE detonation, segment sl-cl in Fig. [H can be viewed as 
a C-detonation overdriven to fc = {Dcj,nse/ Dcj,cY ■ AH C-detonations with velocities 
between Dcj,c and Dcj,nse, such as s2-c2, can be considered as overdriven detonations 
with overdrive parameters 1 < fc < {Dcj,nse/ Dcj^cY ■ Equivalently, they are underdriven 
detonations characterized by / ranging from {Dcj,c / Dcj^nse)"^ to one. The same detonation 
can be characterized by using either / or fc which are related by f|T2|) . In this paper we choose 
to uniformly parametrize all detonations using /. The above consideration equally applies to 
0-detonations with the replacement of qc with qo and Dcj,c with the corresponding value 
of Dcj,o- 

Since C and 0-detonations have lower energy release, their temporal and spatial scales 
are larger than those of the NSE detonations. Fig. [5] shows the reaction zone length xc and 
Xo for C-detonations and 0-detonations as a function of / for densities p = 10^, 3 x 10^, 10^ 
g/cm^. The curves begin on the left at the values of / which correspond to the CJ regimes 
of C and 0-detonations. On the right, all curves terminate at / = 1 where the C and O- 
detonation becomes a part of a complete CJ NSE detonation moving with velocity Dcj^nse- 
Spatial scales of C and 0-detonations increase with decreasing / exponentially. For CJ 
C-detonation xc exceeds the thickness of the / = 1 detonations by a factor ~ 30 — 100. 
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4. One-dimensional propagation of CO detonations 

4.1. Detonation stability 

Numerical simulations of detonation propagation are listed in Table. [TJ Details of a 
numerical approach and test simulations are described in Appendix |X1 For each simulation 
the table gives the initial density po, the overdrive / and other parameters of a steady- 
state detonation. The column seven of the table describes the non-steady behavior of a 
detonation wave found in the simulations: " S" - stable propagation, " O" and " D" - unstable 
propagation. In 0-type cases we observed several cycles of oscillations of the detonation 
wave. In D-type cases the detonation rapidly decayed, the leading shock wave decoupled 
from the reaction zone and the re-ignition of the detonation did not occur over the time of 
numerical integration. This does not mean that the detonation dies and never re-ignites but 
that the re- ignition does not occur on the time scale of at least ~ L/ag, where L is the length 
of the computational domain and is the sound speed. The re-ignition time for D-type 
detonations is estimated in Sect. 14.31 

Fig. in] summarizes the stability properties of the calculated detonations on the po - f 
plane. The boundary of deton ation stabi l ity for NSE detonations in 0.5C+0.5O mixtures has 



been previously calculated in iKhokhlovl ( 1l993l ). It was found that the CJ NSE detonations 



(/ = 1) are stable for po < 2 x 10'' g/cm^. For higher densities the boundary of stability 
passes through overdri ven detona t ion re gimes with / > 1. The stability curve shown in Fig. [6] 



combines the results of IKhokhlovl (119931 ) for pq > 2x 10^ g/ cm^ with the results of this paper 
for Po < 2 X 10^ g/cm^. At low densities the curve passes through overdriven regimes of 
C-detonations with / < 1. The stability curve crosses the / = 1 line at po = 2 x 10^ 
g/cm^. It must be stressed that at all densities the regimes located at the stability curve are 
characterized by /c > 1 and are overdriven with respect to the CJ C-detonation. The freely 
propagating C-detonations are highly unstable and exhibit the D-type behavior. Overdriven 
detonations with / = 1.08 calculated for high densities po > 10^ g/cm^ have the structure of 
the C-burning layer which is very close to that of freely propagating pathological detonations. 
These detonations are highly unstable and exhibit the D-type behavior as well. 

Changes in the detonation behavior when po is fixed and / is decreasing are illustrated in 
Fig. |7]for Po = 3 X 10^ g/cm^. At this density the detonation with / = 0.76 is stable whereas 
the detonation with / = 0.71 is unstable, from which we may deduce that the boundary of 
stability at this density is located between these two values of /. The linear interpolation 
gives fs — 0.73. The detonation with / = 0.71 is mildly unstable (Fig. |7^). The initial slow 
growth of perturbations is followed by regular oscillations with a period 11 ~ 9.5tc and an 
amplitude ~ 0.1 of the steady-state post-shock pressure. For the detonation with / = 0.62 
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the instability develops more rapidly. The quasi-periodic oscillations have a much larger 
period, 11 ~ 50tc, and a much large amplitude ~ 4 of the steady-state post-shock pressure 
(Fig. [7)d). For a detonation with / = 0.52 the initial rapid increase of post-shock pressure is 
followed by a rapid decay of a detonation wave (Fig. [7t). Fig. [TJi compares calculations of the 
/ = 0.55 detonation made with A ~ xc/30 and L = lOOxc (dots) and with A ~ xc/100 and 
L = 640xc (solid line). The comparison shows that the behavior is resolution-independent 
and not influenced by the boundary conditions. 

4.2. Oscillation cycle 

Figs. [HI M and [10] show the distributions of pressure, temperature, and carbon mass 
fraction inside the reaction zone at different moments of time during the third cycle of 
oscillations (Fig. [7)d) for the detonation po = 3 x 10^ g/cm^ and / = 0.62. During the 
minimum of post-shock pressure, times ti to t^, the leading shock is weakened and the post- 
shock temperature is smaller than the corresponding steady-state post-shock temperature. 
As a result the reaction zone is approximately ten times wider than xc- 

An important consequence of the widening of the reaction zone is the decrease of the 
gradients of temperature and reactivity of matter behind the leading shock. Flattening of 
the gradient of reactivity causes the development of a supersonic spontaneous reaction wave 
behind the leading shock between times ^3 and t^. The emergence and amplification of the 
spontaneous wave is clearly visible in Fig. [8] at and ts as a growing peak of pressure 
propagating to the right. Between and the forward part of the spontaneous wave 
steepens into a shock and the spontaneous wave transforms into a secondary detonation 
wave which propagates through the shocked and compressed un-reacted matter towards the 
leading shock. 

Between tg and the secondary detonation overcomes and interacts with the leading 
shock and passes into the cold unburned matter where it becomes the main detonation wave. 
The interaction of the secondary detonation with the leading shock generates a rarefaction 
wave which travels back into the burned material. The rarefaction leads to a gradual weak- 
ening of the leading shock and widening of the reaction zone {ts to ti2)- By the time ti2 the 
leading shock is weakened to a degree that the distribution of physical parameters behind 
the shock resembles the distribution which existed in the beginning of the oscillation cycle 
at time ti. At this moment the new oscillation cycle begins. 



Generation of a spontaneous wave and the subsequent transiti on to a detonation in 
compressed matter occurs via the Zeldovich gradient mechanism (IZeldovich et al.l Il970l : 
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Lee &: MoenI Il980l ) in which the pressure wave is amphfied by the energy release source 
which travels with the spontaneous velocity Dsp = ( '^^^Ji^'*'' ) 5 where t(T) is the energy 



\ dx 

release reaction timescale. For an efficient amplification of a pressure wave, Dgp must be 
comparable to or exceed the local speed of sound, and the spontaneous reaction source must 
act for a sufficient period of time. The diminishing of the temperature gradient caused by 
the weakening of the leading shock favors both these requirements. 

Since r is an exponential function of T, the size of the reaction zone during the minimum 
of the detonation cycle depends exponentially on /. For / = 0.71, the increase of the reaction 
zone at the minimum of oscillations is much less pronounced than for / = 0.62. As a result 
Dgp becomes smaller and the generated pressure wave is weak and does not transition to 
a secondary detonation before it reaches the leading shock. In this case we observe a low 
amplitude, nearly sinusoidal oscillating behavior of a detonation wave. With decreasing / 
the length of the reaction zone at minimum increases dramatically. In the D-case / = 0.55 
the reaction zone rapidly separates from the leading shock and the distance between the 
shock and the reacted matter becomes much larger than the length of the computational 
domain L. In this case we observe the decay of the detonation wave but the re-ignition and 
formation of a spontaneous wave is significantly delayed and is not observed. 



4.3. Re-ignition time 

We now proceed to estimate a characteristic time-scale of re-ignition for D-type deto- 
nations. Fig. [11] illustrates the D-type behavior of a detonation with po = 3 x 10^ g/ cm^ 
and / = 0.55. In this case the initial perturbations grow and then the reaction zone rapidly 
separates from the leading shock. The post-shock pressure decreases from the initial val- 
ues P ^ 13Po (time ti) to P ~ SPq at t-r, and temperature from Tg ~ 2.6 x lO^K to 
Ts ~ 1.8 X lO^K. By the time tr the detonation pressure spike disappears and the pressure 
in the shocked material becomes constant and approximately equal to the pressure Pd of the 
corresponding steady state detonation. 

The re-ignition time can thus be estimated as follows. For a given po and / we first 
calculate the detonation pressure Pd assuming qd = qc- Next, we calculate the post-shock 
temperature T, and post-shock density ps of a shock with the post-shock pressure = Pd- 
Finally, using the values of T,, and the initial carbon composition Yqa we calculate the 
timescale TciTg, p^, Ycfl) for the leading C^^ + C^^ reaction. This time serves as an estimate 
of the reignition time. 

We calculate tc as follows. The kinetic of the C^^ + reaction may be approximated 
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as 



(13) 



where Q ~ 84.165 x lO^K^/^ T„ = T/(l + 0.067 x IQ-^T) and v4(T) is a known weak 
(non-exponential) function of temperature (IFowler et al.lll975[ ). 

In what follows the subscript i indicates initial physical values at the beginning of the 
reaction. The e-folding time of C^^ + C^^ reaction is 



r, 



C,i 



1 dYc,^ 
Yc,i dt 



Rc 
Yc 



(14) 



if we assume that the reaction occurs at constant T = Ti and p = pi. To evaluate the actual 
time-scale tc we must take into account self-acceleration of the C^^ + C^^ reaction caused 
by increase of T. 



The evolution of T is governed by 



dT 
'dt 



Yr 



C,i 



RciT,p,Yc 



c 



(15) 



where qc is the energy release in -|- reaction per uni t mass of CO mixture with initial 
carbon composition Yc^i- Using iFrank-Kamenetskiil (119671 1 (FK) approximation we assume 
that all values in ( IT^ can be evaluated at constant initial conditions except the exponent in 
Rc{T). In FK approximation we can re- write ( IT5l) as 



dT 
'dt 



T- 



exp 



Q Q 



1/3 



T 



where 



T, 



TTd = 



1/3 



exp 



- , T 



(91nT, 



(16) 



(17) 



is the temperature e-folding time evaluated at constant initial conditions and 



6 



3T. 



1/3^ 



Q 



T- 



is the Frank-Kamenetskii factor. Integration of ( !T6|) gives 



Tc - Tt,, 



exp ( e-M ^ - 1 



Ti 



d ( - 



e ■ TTA. 



(19) 



Formally, the upper limit in the integral f lT9|) should be equal to some final temperature 
Tf > Ti reached at the end of burning. Tj may be replaced with oo since the integral in 
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f fT9|) is mostly accumulated at temperatures T c:^ Ti. FK factor O in f|T9l) characterizes 
shortening of the reaction timescale due to self-acceleration of burning. We note that in the 
FK approximation 

so that Tc may also be written as 

Fig. [T2] shows temporal and spatial re-ignition scales as a function of po which were estimated 
for the borderline D-cases (see Table [1] and dashed line in Fig. [6]). Both scales increase 
exponentially with decreasing po- At densities less than ~ 10^ g/cm^ the time scale becomes 
greater than the explosion time of a white dwarf (~ 1 sec) and the spatial scale, estimated 
by multiplying the velocity of the leading shock of the detonation with the reignition time, 
becomes greater than the scale of the exploding white dwarf. At these densities the one- 
dimensional propagation of a C-detonation will be impossible. 



4.4. Stability criterion 

In general, the instability is a consequence of a positive feedback between the hydrody- 
namical perturbations and the temperature-dependent energy release inside the detonation 
wave. The feedback loop consists of three parts. First, the temperature perturbation leads 
to the increase of the reaction rate. This effect depends on the sensitivity of R to T. Second, 
the increase of R leads to the increase of the energy generation rate. This effect is propor- 
tional to the amount of nuclear energy which is released by burning. Third, the increase 
in the energy generation leads to a faster rise of T. The latter effects depends on specific 
heat of matter. For illustration, consider a thermally isolated, uniform fluid element of a 
carbon-oxygen mixture at constant volume. The evolution of T in the element is described 
by (|T^ . To compare the feedback loop for fluid elements with different initial temperatures 
and densities we scale t to the time rc j which is required to burn half of initial carbon at 
constant initial conditions (see ( 11^ ). Then in FK approximation (ITSl) becomes 

^ = r. \ RciT,p,,Yc,) 

where i = t/Tc,i- 

Consider now two identical fluid elements, one with initial temperature Tj and another 
one with perturbed initial temperature Tj + A T, AT << Tj. The evolution of AT with t 
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is given by 

rfAT d{T + AT) dT fdlnTA f qc\ fdlnRi 



dt dt dt V d\n ti J \ ei j \ d\n Ti 

Hence AT(t) oc exp(/3t) where the increment 



AT. (23) 



We see that the increment is the product of three dimensionless factors. First factor in ([2] 
describes sensitivity of R to temperature, second factor characterizes the relative energy 
content of the fuel, and third takes into account specific heat of matter. 

Table [T] gives /3 and individual multipliers in fl2^ for computed detonations. Inspection 
of the table shows that variations of /3 with po and / are mostly caused by variations of 
the second and third multipliers in (!24|) . For example, for p = 10^ g/cm^ the value of (3 
varies by a factor of ~ 2.90 for / varying from 1 to 0.61. qc/^s and {din T/d\ne) vary 
by a factor ~ 1.6 whereas {din R/d\n T) = 0^^ varies by a factor of only ~ 1.1. The 
situation is similar for other densities, as well. The reason for large variations of the second 
and the third multipliers is related to the dependence of internal energy per unit mass and 
specific heat of degenerate matter on p. On the other hand, is density-independent 
(screening corrections at post-shock temperatures are small and can be neglected) and is a 
weak function of temperature, O^^ ~ T^^^^. As a result, relative variations of are small 
and this factor does not contribute to the variations of /3 significantly. 

We note that the relative importance of the factors controlling the instability of a 
detonation is different in C-detonations in supernovae and in terre strial detonations in re- 



active gases. The equation of state for terrestrial gases (for example iBurcat fc Ruscidl2005[ ) 



is such that e,; is usu ally independent of p and increases with T linearly or faster (e.g. 



Burcat fc Ruscidl2005[ ). Specific heat in reactive gases is also density- independent and may 
be a constant or a weak function of T. On the other hand, Arrhenius type reaction rates of 
chemical reactions are R ~ exp(— Q/T) and thus = T/Q is a stronger function of T 
compared to for C-detonations. As a result, the sensitivity of chemical reaction rates 
on T play a significant role for terrestrial detonations. 

Close inspection of Tabled] shows that the transition from stable to unstable behavior of 
a C-detonation with decreasing / wave typically occurs for /9 ~ 9 — 10, virtually independent 
of Po- The borderline between 0-type and D-type detonations corresponds to /3 ~ 14 — 15. 
Of course, the correlation is not perfect. Perfect correlation should not be expected since the 
detonation wave is a distributed system in which physical conditions vary inside the reaction 
zone. Communication of information from one point to another involves interaction of many 
fluid elements and occurs with a finite speed. Derivation of is done within a simple 
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one-zone model and does not take any of this into account. Nevertheless, it seems that the 
derivation catches a significant part of the feedback. 

Finally, the feedback analysis helps to understand why detonation in degenerate carbon- 
oxygen tends to become more unstable at higher densities. The leading part of a complete 
CO detonation can be viewed as a C-detonation overdriven to / > 1. The feedback loop 
consideration applies to this detonation as well. With increasing po specific heat of degenerate 
matter decreases drastically and as a result {din T/dlne) in fl24l) increases by a factor of 
almost ~ 10 (see Tabled]). This leads to a noticeable increase of the strength of the feedback 
loop. As a result, greater overdrive is needed for stabilizing the detonation and the stability 
boundary passes through the overdriven detonation regimes (Sect. 14. II and Fig. [6]). 

5. Discussion and Conclusions 

In this paper we studied basic properties and one- dimensional stability of incomplete 
carbon-oxygen detonations in low density environment of SN la. Calculations of the steady- 
state structure of C-detonations show that a Chapman- Jouguet C-detonation has a velocity 
Dc,cj which is a factor of ~ 20 — 25% smaller than the CJ velocity of a complete NSE 
detonation Dcj,nse- This decrease translates into the decrease of the post-shock temperature 
by a factor ~ 0.7 (Tabled]), and the corresponding increase of the half-reaction thickness 
Xc by a factor ~ 30 — 100 (Fig. |5]). The half- reaction time scale tc — Xc/D increases by 
approximately the same factor. While the increase is appreciable, at all relevant densities 
Xc still remains relatively small compared to the the characteristic scale of a WD, ~ 10^ cm. 

Decreased energy release in C-detonations has a major effect on the detonation stability. 
In contrast to complete detonations which are one-dimensionally stable at low densities, 
C-detonations are highly one-dimensionally unstable. The boundary of stability for C- 
detonations (Fig. |6]) passes through detonation regimes which are overdriven with respect to 
CJ C-detonation. Near the boundary of stability unstable detonations exhibit low-amplitude 
pulsations and have a continuous distribution of physical parameters inside the reaction zone. 
With decreasing overdrive the detonation becomes highly unstable and begins to propagate 
by periodically generating a spontaneous reaction wave inside the primary reaction zone. The 
spontaneous wave transitions to a secondary detonation which overcomes and accelerates the 
leading detonation shock. The subsequent decay and the separation of the leading shock from 
the reaction zone gives rise to a new spontaneous wave and starts the next oscillation cycle. 

Transition from stable to unstable behavior, the level of instability, and the length of the 
detonation cycle depend mainly on specific heat and internal energy of degenerate matter 
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behind the leading shock. The sensitivity of reaction rates to post-shock temperature is less 
important. This is in contrast with idealized Arrhenius detonations whose sta bility is usually 



highl y dependent on the sensitivity of the reaction rate to temperature (e.g. iLee &: Stewart 



19901 ). Stability and detonation cycle of a C-detonation in SN la cannot be understood 



without taking effects of the equation of state of degenerate matter into account. 

Thickness of a steady-state one-dimensional detonation wave (Fig. |5]) can be used to 
estimate the minimal resolution required for a correct description of the detonation structure 
and detonation stability. The absolute minimum of numerical resolution can be no less than 
n ~ 10 — 30 computational cells per xc- The minimal spatial resolution Ar depends on 
background conditions in unburned matter and can be estimated for a given density and 
overdrive by taking xc from Fig. O and dividing it by n. Even for the most favorable case 
of a C J C-detonation at po = 10^ g/ cm^ (the upper left end of the solid line for Ig po = 6 
in Fig. [5]) and n = 10 the required resolution 5 r ^ 10^ cm is significantly less than the 
resolution of current global three-dimensional simulations of SN la. For one-dimensional 
lagrangian simulations we can estimate the relative mass resolution 5q = 6M/Mwd as 

^ ^^PoR'xg ^ ( R \ ' fM^\ / 10 \ / Po \ 

^ nMwD XlQ^cm) \Mwd ) \n ) \lQ^g / cm^ ) ' ^ ' 

where M is the lagrangian mass coordinate and R is the radius of a spherical layer with 
density po at the moment of time when the detonation is passing through the layer. For 
Po = 10^ g/cm^, i? ~ 5 X 10^ cm (see below) and n = 10 we obtain 5q ~ 10~^. 

The period and spatial scale of the detonation cycle of a C-detonation increase with 
decreasing density. For p < 10^ g/cm^ the period becomes larger than the explosion time 
scale, ~ 1 sec, of a white dwarf, which suggests that the influence of the instability on 
propagation and quenching of a C-detonation wave could be significant at these densities. 
This may alter the composition of the outer layers. In particular, quenching at higher 
densities may leave more unburned carbon. From the observational point of view it would 
be important to estimate the velocity range of a SN la envelope which may be affected. This 
depends, of course, on a particular ex plosion model . For illustration, we give estimates for 



two delayed detonation (DD) models of iHofiich et al.l (|2002| ). For a normal supernova we take 



the DD model SP0022.20 with total mass Mwd = 1.346M0 and Ni^^ mass Mm = OMMq. 
In this model the detonation wave, which propagated outward through the WD already pre- 
expanded by deflagration, passed through the layer with density po = 10^ g/cm^ when the 
layer was located at the mass coordinate M = I.SOOMq and radius i? ^ 5 x 10^ cm. At 
this moment of time the layers with density po < 5 x 10^ g/cm^ had radii R > 5.6 x 10^ 
cm. From this we estimate the characteristic scale of density variation at these densities 
as L ~ (^^) — 5 X 10^ cm. The corresponding time-scale of density variation can be 
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estimated by dividing L by the difference in the expansion velocities of the layers ~ 10^ 
cm/sec which gives the characteristic time-scale ~ 0.5 sec. The numbers are consistent with 
the values of 10^ cm and 1 sec used in our general analysis in Sections S] and [31 and in 
this section above. After the explosion and acceleration to free expansion, the layer with 
Po = 10^ g/cm^ was accelerated to ~ 23,000 km at infinity. From this we conclude that 
the one-dimensional instability in a normal delayed detonation explosion may affect the outer 
O.O4M0 or 3% of a SN la envelope with velocities V > 23, 000 km/sec at infinity. For a sub- 
luminous DD model 5P0Z22.8 similar analysis shows that the instability may affect outer 
O.O8M0 or 6% of the envelope with velocities V > 14, 000 km/sec. Thus, the instability may 
affect the spectra and the light curve of these supernovae before and near maximum light. 
We note that the mass resolution of the models was 6q ~ 10~^, i.e., two orders of magnitude 
greater that the resolution required for detecting the detonation instability. 

In conclusion we stress that this study is one-dimensional and we caution that the above 
estimates are preliminary. In addition to a str ong one-dimensional instability the detonation 



will be affected by trarisversa l instability (e.g. iGamezo et al.lll999l ). curvature of the detona- 



tion front (e.g. ISharpdl200ll ). and spatial variations of the background density. Transverse 
instability is caused by the same feedback between the energy generation and temperature 
fluctuations which is responsible for the one-dimensional longitudinal instability. Other ef- 
fects depend on the strength of the feedback, as well. Thus, incomplete C-detonations should 
be expected to be more unstable and more sensitive to the above-mentioned effects than the 
complete detonations. Investigation of the multi-dimensional behavior of C-detonations will 
be reported in the next paper (in preparation). 
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Table 1: Properties of NSE and C-detonations in 0.5C^^ + 0.50^^ mixtures. 
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is in g/cm'^; Tg is in lO^K; D is in lO^cm/s. 
''Detonations with / > 1 are calculated assuming NSE in products of burning (Sect. [X^ . Detonations with 
/ < 1 are C-detonations f Sect. 153)) . 
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A. Numerical method 



A.l. Numerical approach 

The equations ([I]) - (jl]) are integra ted using an ada ptive mesh refinement (AMR) re- 
active fiow fiuid dynamic code ALLA (iKhokhlovl Il998l ). Reaction terms are coupled to 
fiuid dynamics by time-step splitting. During the hydrodynamical sub-step the Euler equa- 
tions are integrated with R = q = using an explicit, directiona l-split, second-order accu- 
rate, Godunov-type conservative sch eme with a Rie mann solver (jColella &: Glaz 1985 ) and 
a monotone VanLeer reconstruction (Ivan Leeiill979l ). Hydrodynamical time step is selected 
using the Courant number = 0.9. During the reaction sub-step the kinetic equations together 
with the equation of energy conservation, differential equations 



— - R — - • 

dt ' dt ^' 



(Al) 



are integrated in each cell using a stiff solver with adjustable sub-cycling to keep the accuracy 
of integration of energy and composition to better than 10~^. 

AMR is carried out at the level of individual cells. Size of computational cells at 
refinement level Imin < ^ < ^max is A(/) = L ■ 2^'; where Imax and Imin are maximum and 
minimum levels of refinement, respectively. During the simulation the mesh is refined around 
the shocks, contact discontinuities, and in regions of large gradients of density, pressure, and 
mole fractions of H^, C^^, O^^, Ae^°, Mg^^ and Si"^^. Detected shocks and discontinuities are 
always refined to maximum level of resolution. Refinement based on chemical composition is 
used for species with concentrati ons Yi > 10~^. Further details of the refinement procedure 
can be found in IKhokhlovl (Il998l ) . 



Calculations are performed in a reference frame moving with a constant velocity of a 
steady-state detonation D in order to keep the detonation on the mesh as long as possible. 
A steady state solution corresponding to a detonation velocity D is mapped onto a compu- 
tational domain of size L. A constant supersonic infiow with u = —D, p = po, P = Pq, and 
Y = Yq is imposed on the right boundary x = L. Zero-gradient boundary conditions are 
imposed at the left boundary x = 0. Calculations were performed with numerical resolu- 
tion of ~ 30 cells per half- reaction length, A ~ xc/30, and L ~ lOOx^ and then repeated 
with two times lower numerical resolution and, when necessary, with up to four times higher 
numerical resolution and L ~ (300 — 1000)a;c' in order to confirm the observed detonation 
behavior. 
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A. 2. Test simulations 



Linear stability and one-dimensional propagation of idealized detonations with the one- 
step Arrhenius kinetics and 7-constant equation of state, 



(A2) 



dt ' dt ^dt' p' 

where Q is the activation energy, q is the total energy release, and 7 = 1.2 have been exten- 
sively studied in the p ast diTpenbeckH lQe^: llee fc Stewa"rtlll99ol : Ifle fc LeelEoOsI : IWilliams et al, 
19961 : ISharpe fc Falldl200d and references therein). We carried out a number of test sim- 
ulations in order to verify that the code can reproduce the stability properties of these 
detonations correctly. 

Fig. [T3h shows pressure Ps as a function of t for an overdriven detonation Q = q = 50 
and / = 1.6. For this case the linear analysis predicts a single unstable mode oc e'^* with the 
increment a = ctr -|- iaj = 0-112 =b iO.789; the corres ponding period of oscillations is 11 = 
27r/cr/ = 7.963 (lErpenbecklll964j : iLee &: Stewartlll990l ). The simulations show initial growth 



of perturbations with ctr ~ 0.1. The amplitude of perturbations saturates due to non- linear 
effects at the level of ^ 0.8 of the steady-state post-shock pressure. The numerical period 
of oscillations 11 ~ 8.42 is larger than the li near p rediction. This is c onsistent with previou s 
numerical simu l ations of iFickett fc David (119791 ) ; iHe fc Led (119951 ) ; IWilliams et al.l (119961 ) ; 
Sharpe fc Falld (120001 ) . The increase of 11 can be attributed to the non-linearity of oscillations 
which manifests itself in a large amplitude of pulsations and in a pronounced asymmetry in 
the shape of Ps{t). According to the linear analysis the transition fr om stable to unstabl e 
behavior of the detonations with Q = q = 50 must occur at / = 1.731 (jLee &: Stewartlll990l ). 
In numerical simulations we find the transition at / ~ 1.75 in a reasonable 1% agreement 
with the linear analysis. 

Figs. [T3r b-d) show results of simulations for CJ detonations with / = 1, g = 50 and 
varying Q. The linear stability analysis predicts the boundary of stability at Q — 25.26 
jHe fc Leel fl995[ ). For Q = 25 the authors predict the stable root a = -0.012 ± zO.494. 
Our simulations show a stable propagation for a detonation with Q = 24 (Fig. [T3b). For 
Q = 25 the simulations (Fig. [T3b ) show quasi-periodic slowly decaying oscillations with 
a ^ — 0.011 -|-i O. 49 in good agre ement with the linear predictions. For Q = 27, interpolation 
in Table II of iHe fc Led (119951 ) gives aj ~ 0.470. The simulations Fig. IT^ gi ve similar 
gj 0.47 in good agreement with the linear analysis and non-linear calculations of iHe fc Lee 
(119951 ) (their Figur e 10b). The instability of a detonation tends to increase w i th dec reasing /. 
In agreement with iHe &: Led (Il995l ) ; IWilliams et al.l (119961 ) ; ISharpe fc Falld (I2OOOI ) we found 
that the propagation of the above detonations becomes highly irregular when / approaches 
/ = 1- 
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Fig. 1. — Schematic representation of NSE and incomplete detonation waves on a pressure 
P - specific volume V = 1/p diagram. O - initial state of matter ahead of the shock. S - 
shock adiabat, NSE - detonation adiabat, C - partial detonation adiabat corresponding to C- 
burning with reduced energy release Qd = Qc < Qnse- si - CJnse - Chapman-Jouguet NSE 
detonation; si - post-shock conditions for the CJ NSE detonation; cl - conditions at the end 
of the C-burning layer of CJ NSE detonation. s3 - CJ^ - Chapman-Jouguet C-detonation, 
s3 - post-shock conditions for the CJ C-detonation. s2 - c2 - a C-detonation overdriven with 
respect to the CJ C-detonation and underdriven with respect to the CJ NSE detonation. 
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Fig. 2. — A Chapman- Jouguct NSE detonation structure for 0.5C+0.5O mixture at po = 3x 
10^ g/cm^ and / = 1. Top - pressure, middle - temperature in lO^K, bottom - mass fractions 
of nuclei (solid curves) and nuclear energy release (dashed line). Distance is normalized to 
a half-reaction length of a C-burning layer. Shock wave is on the left (at — oo in logarithmic 
coordinates) . 
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Fig. 3. — An overdriven detonation structure for 0.5C + 0.5O mixture at po — 10^ g/cm^ 
and / = 1.08. Top - pressure, middle - temperature in lO^K, bottom - mass fractions of 
nuclei (solid curves) and nuclear energy release (dashed line). Distance is normalized to a 
half-reaction length of a C-burning layer. Shock wave is on the left (at — oo in logarithmic 
coordinates) . 
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Fig. 4. — Thickness of reaction layers as a function of density for unsupported NSE det- 
onations in 0.5C + 0.5O mixture, xq - half-reaction thickness of C-burning layer, xq - 
half-reaction thickness of 0-burning layer, xsi - half-reaction thickness of 5'i-burning layer, 
xnse - total thickness of the detonation wave. 
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Fig. 5. — Thickness of C-dctonation {xc, solid lines) and 0-detonation {xq, dashed lines) 
as a function of overdrive / < 1 for densities lg{po) = 6.0, 6.5, 7.0; po is in g/cm^. 
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Fig. 6. — Summary of calculations of time-dependent propagation in 0.5C + 0.5O mixtures. 

Open circles - stable detonations; Open triangles -unstable oscillating detonations; solid 
triangles - unstable detonations for which the oscillation cycle was not calculated. Solid 
line - the estimated stability curve. Dashed line - time-scale of the oscillation cycle exceeds 
~ 200tc where tc is the half-reaction timescale of a steady-state C-detonation. Dash-dotted 
line - location of Chapman- Jouguet detonations, / = 1. 
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Fig. 7. — Post-shock pressure Ps{t) for C-detonations with different overdrives at po = 3x10^ 
g/cm^. (a) / = 0.71; (b) / = 0.62, and (c) / = 0.55. (d) is the early portion of (c) 
and compares two simulations with different distance between the shock and the outflow 
boimdary (dots - shorter distance; see text). Time is normalized to a half- reaction time- 
scale of a steady-state detonation. 



- 31 - 




130 



140 



150 



160 



X/Xr 



Fig. 8. — Oscillation cycle of the / = 0.62 C-detonation at po = 3 x 10^ g/cm^ shown in 
Fig. [Th- Plots of pressure, P{x), for various moments of time during the third cycle. Distance 
is normalized to a half-reaction length-scale of a steady-state detonation. Times are (1) - 
113.2, (2) - 124.3, (3) - 134.1, (4) - 136.2, (5) - 137.5, (6) - 137.9, (7) - 138.8, (8) - 139.2, (9) 
- 139.6, (10) - 142.2, (11) - 174.5, (12) - 183.1; normalized to a half-reaction time-scale of a 
steady-state detonation, (a) - formation of the spontaneous reaction wave and the secondary 
detonation in post-shock matter behind the weakening leading shock (times 1 to 6); (b) - 
re-establishment and subsequent weakening of the primary detonation (times 7 - 12). 




Fig. 9. — Oscillation cycle of the / = 0.62 C-detonation at po = 3 x 10^ g/cm^. Same as 
Fig. Obut shows plots of temperature T{x). 
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Fig. 11. — Propagation of a detonation wave for po = 3 x 10^ g/cm^ and / = 0.55. The 
figure shows the development of the instabihty and separation of the reaction zone from 
the leading shock. Distance is normalized to a half-reaction length-scale of a steady-state 
detonation. Times are (1) - 115.6, (2) - 134.1, (3) - 136.2, (4) - 137.5, (5) - 137.9, (6) - 138.4, 
(7) - 138.8; normalized to a half-reaction time-scale of a steady-state detonation. 




Fig. 12. — Solid circles - estimated reignition length-scale; Open circles - estimated reignition 
time-scale (see Sect. I4.3| l. Horizontal dashed line - characteristic spatial scale of an exploding 
Chandrasekhar-mass CO white dwarf. 
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Fig. 13. — Post-shock pressure as a function of time for Arrhenius detonations, (a) - an 
overdriven detonation with Q = q = 50 and / = 1.6. The last three frames show results 
for a CJ detonation with f = 1, q = 50, and varying activation energy, (h) - Q = 24, (c) - 
Q = 25, and (d) - Q = 27. Time is normalized to a half- reaction time-scale of a steady-state 
detonation. 



